SPRAY COMBUSTION STABILITY PROJECT 


NS6“ 70070 


Dr. San-Mou Jeng”, Principal Investigator 
Mr. Ron J. Litchford^, Co-Investigator 

Center for Advanced Space Propulsion 
UTSI Research Park 
Tullahoma, Tennessee 37388 


ABSTRACT 

This report summarizes research activity on the 
Spray Combustion Stability Project, characterizes ac- 
complishments and current status, and discusses pro- 
jected future work. The purpose is to provide a con- 
cise conceptual overview of the research effort so the 
reader can quickly assimilate the gist of the research 
results and place them within the context of their 
potential impact on liquid rocket engine design tech- 
nology. Therefore, this report does not elaborate on 
many of the detailed technical aspects of the research 
program. Technical issues are addressed in substan- 
tial detail in recent AIAA papers which were pre- 
sented at the 2 6 th Joint Propulsion Conference in July 
1990 and are referenced herein. 

I, Project Synopsis 

The Spray Combustion Stability project is a re- 
search program which is cost shared by CASP and 
Rocketdyne Div., Rockwell Int. Corp. The prin- 
cipal goal of the project is to advance liquid rocket 
engine design technology, specifically analysis of res- 
onant combustion instability. Project work began in 
fiscal year 1989 and is scheduled to continue through 
fiscal year 1991. 

The specific objectives of the research effort as 
elaborated on in the Statement of Work (SOW) is 
to evaluate spray combustion physical submodels for 
application in CFD instability analyses. Submodels 
found deficient were to be refined or replaced by newly 
developed ones. In either case, the resulting submod- 
els were to be validated with experimental data inso- 
far as possible. 

In the original SOW, Rocketdyne ’s ARICC spray 
combustion code was to serve as the testbed for sub- 
model evaluation with improvements incorporated as 
second-generation code developments. During the ini- 
tial part of the first year’s effort, it became apparent 
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that evaluation using the ARICC code was a tedious 
and lengthy procedure. Since a large amount of com- 
puter time and some degree of expertise with the code 
was required to obtain a converged solution, produc- 
tivity was low. A decision was therefore made to eval- 
uate and develop submodels independent of and ex- 
ternal to ARICC. This method proved not only more 
efficient but more practical since physical processes 
could be studied independently. Final implementa- 
tion of submodel improvements is being left to the dis- 
cretion of code developers and users at Rocketdyne. 

Based on a thorough review of the scientific liter- 
ature and some contemplative reasoning, we have thus 
far concentrated our attention on two spray combus- 
tion physical processes which are inherently impor- 
tant in the proper assessment of resonant combus- 
tion instability and which are either inadequately or 
inefficiently simulated by current mathematical sub- 
models. The physical processes considered are (a) 
high-pressure droplet gasification and (b) dispersion 
of the spray droplets and the associated spatial dis- 
tribution of mass and energy release. These processes 
have been the focus for submodel development during 
this project. The technical progress achieved thus far 
is summarized in the following sections. 

II. Technical Discussion 

The guiding rationale in selecting spray combus- 
tion physical submodels for evaluation was the known 
or potential role of a particular physical process as a 
combustion instability mechanism or the pressing re- 
quirement to simulate a particular physical process 
with greater accuracy and efficiency during unstable 
combustion than is normally required during stable 
combustion. The basis of selection led to consider- 
ation of the two physical processes of high-pressure 
gasification and droplet dispersion. Other physical 
processes are certainly fundamental to combustion in- 
stability, but these seemed to have the greatest poten- 
tial for significant improvement within the scope and 
limits of this project. 

Because liquid rocket engines typically operate 
at high chamber pressures, generally above the crit- 
ical pressure of the propellants, the process of high- 
pressure droplet gasification plays a fundamental role 
in combustor dynamics. This role can become even 
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more central when considering acoustic pressure fluc- 
tuations associated with resonant instability. In in- 
stability analyses, it has become a common assump- 
tion to consider propellant vaporization as a rate- 
controlling feedback mechanism for Rayleigh-type 
driving of the acoustic waves. Thus, because the va- 
porization process is demonstrably sensitive to pres- 
sure, proper assessment of droplet gasification can 
be essential to the validity of the instability analy- 
sis. In many cases, this submodeling issue has been 
averted by the application of low-pressure models 
under high-pressure conditions with the hope that 
the low-pressure model sufficiently approximates the 
high-pressure process. Evaluation of high-pressure 
gasification submodeling was therefore seen to be and 
became a major concern of this research program. 

Accurate description of droplet dispersion and 
the associated spatial distribution of mass and en- 
ergy release was also found to be an important phys- 
ical spray process that needed submodeling upgrades 
in obtaining an efficient technique suitable for insta- 
bility analyses. The currently accepted modeling phi- 
losophy is to represent many physical droplets by a 
computational parcel which interacts with the flow as 
a single effective droplet. That is as a point source 
for interphase transport. To simulate dispersion due 
to turbulent fluctuations, a stochastic sampling pro- 
cedure is conventionally utilized. However, the dif- 
ficulty with this approach, particularly for instabil- 
ity analysis, is essentially one of computational effi- 
ciency.. To obtain reasonable dispersion results using 
the current stochastic procedure, a statistically signif- 
icant number of parcels (typically several thousand) 
must be evaluated. Similarly, to prevent the genera- 
tion of large fluctuating interphase transport sources 
in space and time (computational shot noise), a large 
number of parcels must be evaluated to more uni- 
formly distribute the source effects throughout the 
computational cells. The level of accuracy desired 
for instability analyses therefore imposes severe com- 
putational penalties by requiring an excessive num- 
ber of parcels using current stochastic sampling tech- 
niques. Improved modeling of the dispersion process 
was needed and this problem has been addressed in 
our research effort. 

High-Pressure Droplet Gasification 

The high-pressure droplet gasification modeling 
effort has been ongoing since inception of this project. 
Early submodel development has been summarized 
in a 1989 interim report 1 with later developments 
included in a recent AIAA conference paper. 2 The 
reader is referred to these supplementary publications 
for mathematical details and graphical presentation of 
computational results. The intent of this section is 


TABLE 1 

Gasification Modeling Accomplishments 

(1) Formulation of transport model for subcritical 
vaporization including surface regression, bulk 
thermal expansion, and internal circulation ef- 
fects. 

(2) Implementation of a high-pressure thermody- 
namic vapor-liquid equilibria model to account 
for real-gas effects and ambient-gas absorption. 

(3) Evaluation of submodel with experimental data 
for n-heptane droplets vaporizing in hot, high- 
pressure nitrogen gas. 

(4) Calibration of thermodynamic model for 0 2 / H 2 
phase equilibria. 

(5) Analysis of LOX droplets vaporizing in high- 
pressure, hydrogen-rich gas. 


to present a brief conceptual description of the formu- 
lated submodel and a summary of performance evalu- 
ation including current capabilitites and limits. Tasks 
accomplished in the high-pressure droplet gasification 
sub modeling effort are summarized in Table 1. 

The aim of the theoretical development was to 
obtain a realistic physical submodel with reasonable 
computational efficiency. A detailed comprehensive 
analysis is not our purpose. In particular, the de- 
sire was to bring high-pressure effects on thermo- 
dynamic and transport processes into consideration 
while maintaining mathematical simplicity. This re- 
quired reformulation of the basic transport model and 
inclusion of a high-pressure multicomponent thermo- 
dynamic vapor-liquid equilibria model. In addition, 
supercritical gasification had to be considered since 
sufficiently high-pressures can drive the droplet sur- 
face to the thermodynamic critical state. 

The reformulated subcritical transport equa- 
tions for the gas-phase are based on the quasi- 
steady assumption with inclusion of surface regres- 
sion and bulk thermal expansion effects which are ne- 
glected in low-pressure models. Transport is modeled 
as a spherically-symmetric, diffusion-controlled pro- 
cess with empirical corrections for convective effects. 
Liquid-phase energy transport is currently handeled 
using 3 alternative models: (1) infinite conductivity 
model, (2) pure conduction model, and (3) effective 
conductivity model. The effective conductivity model 
is the most sophisticated approach and attempts to 
account for internal circulation effects. 

Because real-gas effects and ambient-gas ab- 
sorption becomes important at higher pressures, a 
high-pressure thermodynamic model has been imple- 
mented to provide the species concentration bound- 
ary condition at the drop surface and to obtain a 
corrected enthalpy of vaporization which can devi- 
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ate substantially for the pure component latent heat 
of vaporization. The model is currently applicable 
to weakly polar molecules and mixtures containing 
quantum gases (eg. hydrogen). There are some exten- 
sions of the model that can handle aqueous mixtures 
which try to account for the large dipole moment of 
the H 2 O molecule, but they are highly approximate. 

Once the droplet surface reaches the critical 
state, a supercritical gasification model is needed. 
Unfortunately, there is as yet no definitive experi- 
mental description of the critical transition and su- 
percritical gasification processes. It is therefore nec- 
essary to resort to ad hoc schemes based on plausible 
presuppositions. The simplest approach is to assume 
flash supercritical transition of the entire droplet, but 
this seems unrealistic since the droplet should actu- 
ally have internal temperature and concentration gra- 
dients. To account for a presumed finite-rate process, 
two alternative schemes are utilized. In one scheme, 
the basic idea is to transform to a constant temper- 
ature boundary condition corresponding to the criti- 
cal mixing point and apply the subcritical diffusion- 
controlled transport relationships. This is considered 
as a representative lower limit for the gasification rate. 
An alternative scheme is to model the process using 
an empirical stripping mode breakup expression. This 
is considered as a representative upper limit for the 
gasification rate. 

A FORTRAN computer code was developed to 
evaluate the proposed submodel. Initially, validation 
calculations were made for comparison with experi- 
mental data on n-heptane droplets vaporizing in high- 
pressure nitrogen gas. Agreement was good verify- 
ing the basic validity of the code and the subcritical 
transport model. Because of unavailable experimen- 
tal data, there was no way to evaluate the supercrit- 
ical gasification schemes. 

The code was further modified to consider LOX 
droplets vaporizing in high-pressure, hydrogen-rich 
gas since this propellant combination is of great 
interest for high-performance liquid rocket engines. 
The thermodynamic model was calibrated using 
high-pressure solid oxygen/hydrogen equilibria data. 
High-pressure effects on thermodynamic and thermo- 
physical properties were also considered. Represen- 
tative computations were then made to demonstrate 
the transition process from steady-state vaporization 
to fully transient heat-up to the critical state with 
subsequent supercritical vaporization. Presentation 
and discussion of the results are in an accompanying 
AIAA paper. 2 

Additional validation with experimental data 
would be desirable. It should also be stated that for- 
mulations for including an envelope flame have been 
constructed but not included in the code because 


of the additional complications which would result 
from the required extension of the thermodynamic 
model at this time. Extension of the thermodynamic 
model for a wider range of substances including polar 
molecules (eg., H 2 O) will be necessary before prac- 
tical implementation in a comprehensive combustor 
code. Nevertheless, the model is currently useful for 
examining the general effect of pressure on the vapor- 
ization process and the resulting influence on combus- 
tion instability. 

Turbulent Droplet Dispersion 

The turbulent droplet dispersion submodeling ef- 
fort began in the latter part of 1989. Based on the 
difficulties encountered with current dispersion mod- 
els as previously discussed, the conceptual basis for 
a novel modeling technique aimed at avoiding these 
deficiencies was conceived. Development of the model 
has proceeded at a quick pace, and technical progress 
has been summarized in a recent AIAA conference 
paper. 3 The reader is referred to this publication for 
mathematical details and presentation of computa- 
tional results. This section will be devoted to pro- 
viding a conceptual description of the model and a 
report of model evaluation efforts to date. Tasks ac- 
complished in the droplet dispersion submodeling ef- 
fort are summarized in Table 2. 

The fundamental concept of this new model was 
to combine a statistical spatial distribution represen- 
tation for each parcel with a statistical turbulent dis- 
persion width transport scheme. Specifically, each 
parcel represents a group of physical droplets through 
a normal probability density function (PDF) in space. 
The instantaneous mean position of each parcel PDF 
is determined from Lagrangian tracking using either 
a deterministic or stochastic sampling procedure. In 
addition the variance of each parcel PDF is deter- 
mined from a statistical, turbulence-induced mean- 
squared dispersion which depends upon previous par- 
cel interactions with the characteristic turbulence 

TABLE 2 

Dispersion Modeling Accomplishments 

(1) Conceptualization and mathematical formula- 
tion of a novel model for turbulent droplet disper- 
sion with PDF representation for droplet parcels. 

(2) Proof-of-concept computations for new model in 
comparison with conventional point source model 
for nonevaporating particles injected into nearly 
homogeneous turbulent flow. 

(3) Proof-of-concept computations for new model in 
comparison with conventional point source model 
for nonevaporating particles injected into a round 
turbulent jet. 
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properties of the flow. This approach allows for a 
smooth distribution of source terms over the grid cells 
with reduced sampling requirements for obtaining the 
dispersion width. A concept schematic of the model 
with stochastic sampling for the mean parcel position 
is depicted in Figure 1. 
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Figure 1. Concept schematic of proposed dis- 
persion width transport model with stochastic 
sampling for mean parcel position. 

The dispersion width transport scheme is the key 
to increased efficiency since it makes use of statistical 
generalizations allowing reduction in sampling. The 
idea is to retain the concept of parcels interacting with 
a sequence of turbulent eddies as in the stochastic 
point source method but to generalize it by formulat- 
ing the mathematics in such a way to allow statistical 
treatment. This is accomplished by repeated appli- 
cation of the linearized particle equations of motion 
for each eddy interaction with successive substitution 
such that a single expression is obtained describing 
dispersion in terms of the turbulent velocity fluctua- 
tions in each eddy, parcel/eddy interaction times, and 
the inertial time lag for each interaction. By introduc- 


tion of a characteristic turbulent velocity fluctuation, 
it then becomes possible to express the mean-squared 
dispersion as the variance of the PDF after the n th 
eddy in terms of preceeding parcel/eddy interaction 
properties. With the mean and variance of each parcel 
PDF, it is possible to combine the PDF distributions 
to obtain the physical particle distribution profile. 

After following through with the described math- 
ematical formulation procedure, model evaluation 
was considered. Two simple test cases chosen were 
for nonevaporating particles injected into nearly- 
homogeneous grid-generated turbulence and into a 
round, turbulent air jet. In both cases particle load- 
ing was small and interphase momentum coupling was 
neglected. The new model was to be compared with 
the stochastic point source model which has been pre- 
viously validated in the scientific literature. FOR- 
TRAN computer codes were developed for these pur- 
poses and computations carried out. Presentation 
and lengthy discussion of the results are in the ac- 
companying AIAA paper. 3 

In summarizing these results, the new disper- 
sion width transport model was found to have sig- 
nificantly improved computational efficiency over the 
stochastic point source method while producing good 
approximations for the physical particle spatial dis- 
tribution profiles. Two schemes were considered for 
evaluating the mean parcel position: (a) Determinis- 
tic Dispersion Width Transport Model (DDWT) and 
(b) Stochastic Dispersion Width Transport Model 
(SDWT). For nearly-homogeneous turbulence, deter- 
ministic tracking with a single parcel (DDWT model) 
was found to sufficiently sample the turbulence prop- 
erties and produce accurate dispersion results as 
shown in Figure 2. For the jet flow with nonhomo- 
geneous turbulence, however, deterministic tracking 
was inadequate. More thorough sampling of the tur- 
bulence was required through stochastic tracking of 
the mean parcel position (SDWT model). Computa- 
tional efficiency was still impressive, though, as dis- 
persion results using no more than 200 parcels were 
found to produce a good approximation to dispersion 
results using 20,000 particles in the stochastic point 
source methodology. An example of such results is 
shown in Figure 3. 

The potential impact of this new dispersion 
model on the accuracy and efficiency of spray combus- 
tion CFD codes, particularly with regard to instabil- 
ity analyses, is significant. Extension of the model to 
include evaporating droplets with interphase coupling 
of mass, momentum, and energy is obvious follow-up 
work for the future. It appears that additional adap- 
tation for considering acoustic perturbations and im- 
provements in efficiency are also possible, and model 
development is continuing. 
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Figure 2. Proof-of-concept computations for 
the deterministic dispersion width transport 
model. 

III. Project Direction 

Research emphasis during the final phases of this 
project will be on development of the droplet disper- 
sion model. The potential impact of this method on 
spray combustion modeling is sufficiently great to jus- 
tify concentrated effort. Further development of the 
high-pressure droplet gasification model will require 
new experimental data which probably won’t be avail- 
able during the project lifetime. However, limited ef- 
fort will be devoted to assessing the general influence 
of pressure on vaporization-controlled combustion in- 
stability using the current model. 

Tentative plans for systematic development of 
the droplet dispersion model have been made. Work 
has already began on extension of the model to 
include evaporation effects.. Based on experimen- 
tal work in the literature, evaluation of vaporiz- 
ing monodispersed methanol droplets in a methane- 
fuelled jet diffusion flame is in progress. Eventually, 
interphase transport coupling will be addressed in an 
effective manner. Other plans include improving the 
efficiency of the procedure by limiting the number of 
interactions allowed to influence dispersion in a given 
eddy. Extension to recirculating flows and adaptation 
for acoustic perturbations are also under considera- 
tion at this time. 



0.0 0.5 1.0 1.5 2.0 2.5 3 0 

r la 


Figure 3. Proof-of-concept computations 
for the stochastic dispersion width transport 
model. 

A simplified assessment of the high-pressure va- 
porization process on combustion instability is be- 
ing undertaken concurrent with the droplet disper- 
sion model development. The basic concept of this 
approach is to apply a simple analysis of the insta- 
bility response factor using low-pressure and high- 
pressure vaporization models to quantify the effect 
of pressure on the feedback mechanism. Such an 
approach has practical application as a combustion 
response module for combustion instability design 
codes. The droplet vaporization code has been mod- 
ified to simulate a simplified 1-D combustor with 
n-heptane droplets vaporizing in the stoichiometric 
products of combustion with oxygen. Sinusoidal os- 
cillations characteristic of acoustic mode instability 
in a cylindrical chamber are superimposed on the 
combustor flow. The high-pressure thermodynamic 
model takes account of real-gas effects including H->0 
molecules using an approximate engineering model. 
Absorption of ambient-gas into the droplet is ne- 
glected, however. The analysis is limited to subcriti- 
cal vaporization since supercritical gasification mod- 
els are dubious. A total of 36 droplets are injected 
per cycle of oscillation with the resulting vaporization 
histories combined to obtain the instability response 
factors based on the Rayleigh criteria. The response 
factors are to be evaluated for a given mean cham- 
ber pressure as a function of frequency. Comparison 
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of results for low-pressure and high-pressure vapor- 
ization models will yield a measure of importance of 
high-pressure vaporization effects on resonant com- 
bustion instability. The greater the response factor 
the greater the tendency to drive instability. Prelim- 
inary computations at moderate pressures indicate 
that consideration of high-pressure effects implies a 
substantial increase in the response factor. This in- 
crease is expected to be even larger at higher pres- 
sures. 


IV. Concluding Remarks 

This report has attempted to describe the accom- 
plishments and current status of the Spray Combus- 
tion Stability Project as well as indicate the planned 
direction of future research. While research effort in 
this project has been primarily restricted to physi- 
cal submodel development, it is desirable that these 
submodeling concepts be pursued to actual imple- 
mentation into comprehensive CFD spray combustion 
codes. This will hopefully make full scale instabil- 
ity analyses more practicable and useful and improve 
spray combustor design technology. 
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